Absence of topological insulator phases in non-Hermitian PT-symmetric Hamiltonians 
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In this work we consider a generalization of the symmetry classification of topological insulators to 
non-Hermitian Hamiltonians which satisfy a combined PT-symmetry (parity and time-reversal). We 
show via examples, and explicit bulk and boundary state proofs that the typical paradigm of forming 
topological insulator states from Dirac Hamiltonians is not compatible with the construction of non- 
Hermitian PT-symmetric Hamiltonians. The topological insulator states are PT-breaking phases and 
have energy spectra which are complex (not real) and thus are not consistent quantum theories. 
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With recent interest generated in the field of topo- 
logical insulators and superconductorsfT-^, the sym- 
metry classification of (primarily free-fermion) Hamil- 
tonians has resurf aced ll9tiT3l ■ The first example of a 
topological insulator (TI), the integer quantum Hall ef- 
fect (IQHE), is gapped in the bulk and exhibits gapless, 
robust states on its boundaries. Most other examples 
of topological insulatorsfT, 'S'-Sl share these character- 
istic features. One interesting distinction between the 
IQHE and, for example, the quantum spin Hall effect 
(QSHE), is that the IQHE is completely robust to any 
type of Hamiltonian perturbation, while the QSHE is 
only robust to perturbations that preserve time-reversal 
symmetry T. A full symmetry protected classification of 
topological insulators[ll 12J and superconductors [lT| 
based on charge-conjugation, time-reversal, and chiral 
symmetries was then unified into a periodic table llT3]| . 

In addition to the classification theory of TI's, Hamil- 
tonian symmetries are important in the theory of 
non-Hermitian Hamiltonians IITIHTtI where it has been 
shown that non-Hermitian Hamiltonians can still de- 
scribe viable quantum systems as long as PT symme- 
try is unbroken. By PT we mean the combined opera- 
tion of a parity/ inversion symmetry P, and T. Given a 
non-Hermitian Hamiltonian H that satisfies [H, PT] = 
one can provide necessary and sufficient conditions that 
the energy spectrum of H be real] 16 1. It is thus natural 
to attempt to extend the periodic table of TI's to non- 
Hermitian Hamiltonians which satisfy i-T-symmetry 
Although our initial hope was to find non-Hermitian ex- 
amples of TI states we instead show that TI phases are 
incompatible with the PT-symmetric construction of (at 
least a large class of) non-Hermitian Hamiltonians. The 
TI states are PT-breaking and exhibit imaginary eigen- 
values even when H has been constructed to preserve 
PT-symmetry. 

In this Letter we start by showing a few pedagogi- 
cal examples of PT-symmetric Dirac Hamiltonians that 
are perturbed away from being Hermitian. We offer 
these examples to show how the construction of non- 
Hermitian PT-symmetric TI states fails. Our focus is on 
Dirac Hamiltonians because they are the minimal mod- 
els for topological insulators. An extension to generic 



insulator Hamiltonians is straight-forward. After the 
examples, we provide a more generic proof of neces- 
sary conditions for the bulk spectra of such Dirac Hamil- 
tonians to have a fully real spectrum {i.e. eigenvalues 
are real for all values of the momenta). Finally, we 
give some arguments about the properties of the gap- 
less boundary states that show that PT-symmetry, non- 
Hermiticity, and topological insulator states do not seem 
to be compatible. 

Before we begin with the examples let us list the rel- 
evant symmetry properties we will use in this Letter. 
The T-operator is represented by T = UK where U 
is a unitary operator and K is complex conjugation. 
Depending on U we can have T^ = ±1. For a Bloch 
Hamiltonian H{p) to be invariant under T we must 
have TH{p)T^^ = i/(— p). The P-symmetry which we 
will call parity when needed, is a unitary operator with 
P^ = +1. There is no requirement on which spatial co- 
ordinates the P operator inverts and for now we will 
leave it generic. The condition that a Bloch Hamilto- 
nian be P-invariant is PH{p)P^^ = H{p) where p is 
a symbol characterizing a given P by indicating which 
coordinates are inverted and which remain unaffected. 
As an example, if all the coordinates are inverted then 
p = —p. Finally, the condition that a Bloch Hamil- 
tonian be PT-symmetric is PTH(p){PT)-^ = H{-p) 
and equivalently requires that H{p) is either both odd 
or both even under P and T separately. We occasion- 
ally mention charge-conjugation symmetry C which re- 
quires CH{p)C-^ = -H*{-p). 



ID DIRAC HAMILTONIANS 

We will begin the examples in ID with the gapped 
continuum Dirac Hamiltonian 



Hid{p) = vppa^ 



(1) 



which is P, T, and PT symmetric with the symmetry op- 
erators T = K, P = a^ and p = —p. This Hamiltonian 
also satisfies a C symmetry with C = a^ and provides 
an example of a Z2 topological insulator protected by 
C symmetry! 12 1 . We can add an additional term to the 
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FIG. 1. Energy spectrum of H[o,iatuce. (Eq- [5} with oper\ 
boixndaries vs. m (a)A = gives a fully real spectrum with 
zero-modes for —2 < m < 0.(b) real part of energy spectrum 
and (c) imaginary part of energy spectrum for A — 0.1. Com- 
plex eigenvalues exist in the former topological phase. 



Hamiltonian to get 



Hid{v)^P(^^ 



iAa'' 



(2) 



which is P-odd, T-odd, C-even, PT-even and non- 
Hermitian and where we set vp = I- Since we have 
not broken C-symmetry the Z2 classification naively 
should remain intact. ID Dirac Hamiltonians with non- 
Hermitian potentials have been studied in, for example, 
Ref . 18 The energy spectrum of this Hamiltonian is sim- 
ple to calculate: 



s± = ±Vi 



pz 



-A2 



(3) 



which is real as long as jrnj > |A|. So we see it is pos- 
sible to add a non-Hermitian perturbation to the Dirac 
Hamiltonian and keep the spectrum to be entirely real 
for values of m both < and > 0. Note that the per- 
turbation has the special property {iAc^, i/i£)(p)} = 
which will become important later. 

Although it is promising that there is a regime where 
this non-Hermitian Hamiltonian will have a real spec- 
trum, there is already something worrying about the 
spectrum. With A = this model has a gap-closing 
phase transition at ttt. = which separates a trivial insu- 
lator phase from a topological insulator phase|[12|. We 
see here that if A 7^ this phase transition becomes 
destabilized when < |to| < |A|. To properly describe 



the TI phase we need to use a lattice version of this Dirac 
Hamiltonian with a Bloch form 



H[ 



D. lattice 



(p) = (sinp)cr^ + (l + TO-cosp)cr''+iAcr^ (4) 



where we have set the lattice constant a = 1. When 
A = the lattice Hamiltonian is in a trivial insulator 
phase when m < —2 or m > 0, and a topological insula- 
tor when — 2 < m < 0. With non-vanishing A we have 
E± = ±^l + (1 + m)2 - A2 - 2(1 + to) cosp. To be a 
viable spectrum this must be real for all —n < p < tt. 
In Fig. [it we show the spectrum for A = as a func- 
tion of m with open boundary conditions. This clearly 
shows the boundary state zero modes which persist 
when —2 < ttt. < 0. In Fig. [lb,c we show the real and 
imaginary parts of the energy spectrum when A = 0.1. 
This shows that as soon as the system nears the phase 
boundary to the topological insulator state the spectrum 
develops imaginary pieces and thus generates a PT- 
breaking phase. This is a common feature and shows 
that the ID topological insulator state here cannot be 
reached when A 7^ 0. 



2D DIRAC HAMILTONIANS 

2D becomes more complicated because of two rea- 
sons (i) there are two natural definitions of P: parity 
{x, y) — >■ (x, — y) and inversion {x, y) — > {~x, —y) which 
both satisfy P^ = 1 (ii) from the conventional classifica- 
tion theory of topological insulators it is natural to look 
at both 2-band and 4-band models. We won't exhaust all 
these cases, but only provide some instructive examples. 

2-band models: Let us consider the Hamiltonian 
which represents the continuum model of a Chern 
insulator IIT9I : 



^2d(p) = Px<j"' +Pyay + ma^ 



(5) 



with T — K and P = a^. This P corresponds to a par- 
ity symmetry and sends py — > —py so p = {px, —Py)- 
Unfortunately this Hamiltonian cannot be made PT- 
symmetric even before we perturb it because the mass 
term breaks PT. Next, let us consider the same Hamil- 
tonian with T = K and P ~ cr^ where this P is an in- 
version symmetry with p = —p. Again even the base 
Hamiltonian is not PT-symmetric because the Px term 
breaks PT. Finally, if we consider the same Hamiltonian 
with T = ia^K and P — a^ and p = — p the mass term 
breaks PT. From these few attempts it seems like we 
cannot get any interesting non-Hermitian TI Hamilto- 
nians with 2-bands. We will see why this is so in the 
general proof section, but the basic idea is that there 
are no additional matrices M which anticommute with 



4-band models: Let us start with a continuum QSH 
HamiltonianimiSl: 



^2D(p)=p-ri+pyr2 + mr° 



(6) 



r2 



T«, r-^ 



ry (E) I, r^ 



with pi = T^ «) cr^ 
T^ (g) cr^, and r" = T^ ® I. This has symmetry generators 
T = ia^K and P = r° and p = —p. We can perturb 
this Hamiltonian with the non-Hermitian PT-symmetric 
matrices {ir^ (8) I, ir^ cr% ir^ (g) ct\ il (g) cr'} but none 
of these anti-commute with the Hamiltonian and will 
lead to imaginary bulk eigenvalues as we will see in the 
following section. 

We can write down another Hamiltonian in 2D 



^2DB(p)=P:^r2+p^r3 + TOr° 



(7) 



This has symmetry operators T = K and P = F" with 
p = —p. There are two interesting terms with which we 
can perturb the Hamiltonian 



AiJ^t] = iAiPi 



{AaT^ 



(8) 



Adding these terms to the Hamiltonian gives us an en- 
ergy spectrum 



E±=±Jpl+p: 



22 - Af - A2 



(9) 



which is real as long as m > Ag + A4. However, this 
model, with T ^ K is not known to exhibit a robust 
topological insulator state anjrway, but we see there is 
the same problem with the gap-closing transition lead- 
ing to a PT-breaking region. 



NECESSARY CONDITIONS FOR REAL BULK 
EIGENVALUES 

The generic features of the example models is that if 
there is a PT-symmetric, non-Hermitian matrix that an- 
ticommutes with H{p) then we can have real eigenval- 
ues, but the topological insulator phase is PT-breakrng. 
If the only PT-symmetric non-Hermitian matrices com- 
mute with at least one term of H{p) then even the bulk 
eigenvalues will not be real for all of momentum space. 
We formalize this statement now. 

Theorem: For a PT symmetric Hamiltonian of the 
Dirac form H = pT + niQ + iAS where {Q, F} = and 
S, Q, F are matrices with S*^ = Q^ = F^ = 1, the quan- 
tities {Q, S} and {F, S} must vanish to have a real bulk 
spectrum. 

Proof: Assume we have a Hamiltonian 



H{p) =pr + mQ + iAS. 



(10) 



By contradiction let us assume that [5*, Q] = so that 
the anti-Hermitian term commutes with the mass term. 



For the point p = the Hamiltonian reduces to -ff (0) = 
mQ + iAS. Since S and Q commute they can be simul- 
taneously diagonalized. Since S^ ^ Q^ ^ +1 the eigen- 
values of these matrices are ±1. Thus the eigenvalues of 
H{0) can only be ±m ± iA and are always imaginary 
if A 7^ 0. This is a contradiction and so we know that 
{S,Q} = must hold if the spectrum is to be real for all 
allowed p. 

Now we want to prove the second necessary condi- 
tion by contradiction. Assume that {S,Q} = 0, which 
we now know must hold, but [S, F] = 0. The Hamil- 
tonian at finite p is as above and H^ {p) = (p^ + m^ — 
A2)I + ipA{F, 5"}. From our assumptions {F, S} = 2TS 
and (FS*)^ — +1 so the eigenvalues of {F, S"} are ±1. 
Thus the eigenvalues of H^ (p) are 



El = ip' + 
=^ ±E± 



i^ - A^) ± 2ipA 
±v/(p±»A)2 +to2 



(11) 



This means that for all non-zero p the eigenvalues will 
be imaginary. This is a contradiction and thus we have 
proven that both {S, Q} — {S, F} = are necessary con- 
ditions for the bulk spectra to be real. ■ 

Corollary: For generic Hamiltonians of the form 
Hip) = Par^+mQ+iAS* where {F"^, Q} = 0and(F")2 ^ 
g2 ^ 5-2 ::^ 1 we must have {S, T"} = {S, Q} = for all 
values of a. 

Proof: The case involving the mass term is unchanged 
from above. For the momentum term, by contradiction, 
assume that there is a value a = gq such that [S, F""] — 0. 
We can then set all pa = for a ^ uq and use the same 
theorem as above to prove the result. ■ 

Thus, we see that to give a real spectrum we must 
have an S that is a Clifford algebra generator along with 
the generators F"^ and Q. To satisfy PT symmetry we 
need 



{PT)iS{PT)- 



-i{PU)S*U-^p-^ ^iS 
{PU)S* + S{PU) = 



(12) 



this means that if S* — S* then PT symmetry requires 
{PU, S"} = 0. If 5 = -S* then PT symmetry requires 
[PU, S] — 0. The generic energy spectrum for such an 



anti-commuting S will be E± = ±\/J2aPl 



A2. 



BOUNDARY/INTERFACE STATES 

Let us return back to the ID Hamiltonian specified in 
Eqs. 1 2 There is a region of the TI phase {i.e. where 
— 2 < TO < 0) where \m\ could be greater than |A| and 
thus could have a real bulk spectrum, but from Fig. [lb,c 
we see that the phase is still PT-breaking. To illustrate 
why the TI state is PT-breaking we need to consider the 
gapless boundary states. We will do this in the simplest 
possible way and capture the essential details by consid- 
ering two interfaces between regions described by Dirac 



Hamiltonians with m{x) — —ttiq for x < —xq next to 
a region with m{x) — rriQ > for —xq < x < xq and 
finally with m{x) = —mo for x > xq. We assume that 
I mo I > |A| so that the bulk energy spectrum is real in 
each region and that 2x0 ^ hvp/niQ so that for our pur- 
poses the interfaces are independent of each other Note 
that this potential satisfies m{x) = m*{—x) and is thus 
PT symmetric. The ID Hamiltonian is 



H = -i — a^ 
ax 



ra{x)(j'' 



(13) 



This Hamiltonian has two zero-energy boundstates EOl 
one per interface domain wall given by 



*o^(^) 



V2 



exp 



± 



m{x')dx' 



±Xo 



1 

Tl 



(14) 



We see something interesting that occurs here, 

PT^+{x) = P^+*{x) = cT^^+*i-x) = ^^(x) (15) 

This means that ^Pg is not an eigenstate of PT and thus 
we do not expect it to have a real eigenvalue if we add 
a non-Hermitian PT symmetric termllS", "161. The PT 
symmetry transforms one boundary state into the other 
one. This is not so strange since the P symmetry should 
interchange the two ends of the ID system, or in this 
case the states on each domain wall. 

Now let us perturb the Hamiltonian by iAct^, as in 
the ID example above, and focus near a single domain 
wall. The unperturbed system has ^o as a zero energy 
state. Adding the perturbation is trivial since ^o is al- 
ready an eigenstate of a^. Thus we see that perturbing 
the domain wall states changes the energies of the states 
from E — to E ^ ±iA w^hich is obviously imaginary. 
Thus, the domain wall, and by analogy, TI boundary 
states break the PT-symmetry and are not compatible 
with a real energy spectrum. 

To show that this is not a pathological case for the ID 
Hamiltonian let us consider the 2D Hamiltonian in Eq. U\ 
which will still have (unprotected) bound states on mass 
domain walls. We will assume the domain walls are in 
the y-direction this time. The Hamiltonian is then 



H. 



(4) 
2D 



Px 



ay 



3 



TO(y)r" 



(16) 



where p^ is a number. For p^ = we can use the same 
ansatz as in the ID case but this time there are two zero- 
energy solutions per wall. For an upward stepping do- 
main wall we find 



*0A = 



1 

V2 



( 1 \ 



1 



-1 





*0B = 



V J 



1 

72 



/o\ 



(17) 





1 



For Px ^ Q these solutions do not automatically diago 
nalize F^ and we are left with a reduced 2x2 problem 



H\. 



(*o,b=.r2|*oj)=P.<- 



where i,j = A,B. Now if we add the allowed, anti- 
commuting non-Hermitian terms iAiF^ + iA4r"' the ef- 
fective Hamiltonian of the edge states becomes 



J^ij —Px"i 

V 

= Pxcr- 



*od*Airi + zA4r4|*oj) 

zAi<.+*A4< (19) 



This Hamiltonian is simple to diagonalize has energies 



E±=±Jpl-Al-Al 



(20) 



which are imaginary at least at p^; =0 if either of the 
Ai 7^ 0. Thus again the boundary spectrum is imaginary. 
The states ^oa/b are not eigenstates of the PT operator 
and will be transformed into the allowed states on the 
opposite wall by PT. These arguments generalize to the 
other topological insulator classes and show that such a 
phase is generically PT-breaking. 

With this proof in place one can immediately begin 
to search for exceptions. The first place to look would 
be topological superconductors whose boundary states 
have weight on both boundaries. We performed a cur- 
sory test of some classes of topological superconduc- 
tors and were not able to construct an interesting non- 
Hermitian PT-symmetric phase. The presence of Majo- 
rana boundary fermions w^hich are non-local may pro- 
vide a way around the boundary state problem but we 
leave this open for future work. Also, our result does not 
immediately apply to topological insulator states pro- 
tected by point-group symmetries like inversion ll21l423l 
or even C4 symmetry[24] . In this case, however, we be- 
lieve that our result can easily be adapted to at least a 
large class of these models since again they can usually 
expressed using minimal Dirac-type models w^here our 
results would immediately carry over. This does leave 
open the possibility of finding non-Hermitian topologi- 
cal phases, but without Dirac-type Hamiltonians at their 
foundation. 
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